Exponentially decaying correlations in a gas of strongly interacting spin-polarized ID fermions with 

zero-range interactions 
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We consider the single particle correlations and momentum distributions in a gas of strongly interacting 
spinless ID fermions with zero-range interactions. This system represents a fermionic version of the Tonks- 
Girardeau gas of impenetrable bosons as it can be mapped to a system of noninteracting ID bosons. We use 
this duality to show that the T = single particle correlations exhibit an exponential decay with distance. This 
strongly interacting system is experimentally accessible using ultracold atoms and has a Lorentzian momentum 
distribution at large momentum whose width is given by the linear density. 
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One dimensional (ID) many body systems provide a beau- 
tiful example of duality. In this context, duality refers to the 
precise mappings that exist between weakly interacting ID 
bosons and strongly interacting ID fermions and vice versa. A 
number of classic examples of duality in ID systems exist; i) 
the Luttinger liquid 1 1, 2] which is dual to a collection of non- 
interacting bosonic collective modes, ii) the Tonks-Girardeau 
gas of impenetrable bosons |3, 4], which is dual to a nonin- 
teracting gas of fermions and iii) the quantum sine-Gordon 
model, which is dual to the massive Thirring model (5|. 

More recently, a number of groups l"^, 8, 9, 10] have 
introduced equivalent models of spinless ID fermions with 
zero-range interactions. Again, duality appears in these mod- 
els as they can all be mapped to a system of 5-function inter- 
acting ID bosons. These models contain the familiar bosonic 
Tonks-Girardeau gas 01 0] as a specific case as well as a 
fermionic version of the Tonks-Girardeau (TG) gas in 
which strongly interacting fermions are related to noninteract- 
ing bosons. In this paper we discuss the correlations among 
these ID fermions with strong interactions. 

Ultracold atom gases provide an excellent arena in which 
to study ID many body systems with short-range interac- 
tions. Already, ultracold atoms have been confined in quasi- 
ID (highly elongated) harmonic traps fTHl where phase fluc- 
tuations in quasi- ID Bose-Einstein condensates lizlllll and 
the TG gas 1, 14.. .15.1 have been observed. 

Although such experiments are actually three dimensional 
(3D) in nature, their low energy properties can be modeled 
by an effective ID Hamiltonian, as first demonstrated by Ol- 
shanii 1 16]. For quasi-lD bosons, this model consists of mass 
m ID bosons interacting through the potential = ^fo5(z), 
where the interaction strength, g^^, is chosen to correctly re- 
produce the low energy behavior of the full 3D Hamiltonian of 
the harmonically confined atoms. Most importantly, can 
be tuned to any value (even infinite) by changing ratio as/a±, 
where Aj is the 3D i-wave scattering length of the bosons 
and a± is the length scale of the tight harmonic confinement. 
Thus, ultracold atoms in quasi- ID traps can be used to probe 
the well studied uniform ID Bose gas with 5-function interac- 
tions ] 17, 18, 19]. The dimensionless interaction strength for 

this system is Jb = where n is the linear density. In the 



Of critical importance in these studies are the static cor- 
relation functions. These determine experimentally measur- 
able properties such as momentum distributions ] 14] and in- 
elastic collision rates ]20]. For the 5-interacting Bose gas, 
these correlation functions are known in a number of different 
regimes. In the TG regime (Jb — > °°) duality can be used to 
calculate the single particle correlations for all densities and 
distances ]16, 21]. For finite interactions Jb the single- and 
local multi- particle correlation functions are known analyti- 
cally in a number of limits (see i22ll23ll24ll and the references 
therein). Lastly, numerical calculations have covered all inter- 
action strengths for both uniform and harmonically confined 
systems ]25]. The most important feature in all these studies is 
that single particle coiTelations at zero temperature {T ~ 0) ex- 
hibit a power law decay with distance. This power law decay 
of coiTelations is also shared by the classic Luttinger liquid 
model of spinless ID fermions ]26]. 

With this background we turn our attention to the single 
particle correlations in the recently proposed model of spin- 
less ID fermions with zero-range interactions. Our focus is 
upon strongly interacting ID fermions, which are dual to non- 
interacting bosons; a sort of fermionic version of the TG gas. 
We emphasize that this system does not appear to be another 
obvious manifestation of the Luttinger liquid. This becomes 
evident in the nature of the single particle correlations, which 
we find decay exponentially with distance at T = 0. 

Our method consists of using duality to obtain a general 
relationship between the single particle density matrices of 
the fermions and their dual bosons. Thus, once the single 
particle correlations are known for noninteracting bosons in 
any external potential, the fermionic correlations are known 
immediately. We apply these general results to uniform and 
harmonically confined ID gases, and give the experimentally 
relevant momentum distributions for each. 

Spin polarized fermions. Here we introduce and summa- 
rize the main features of the ID spinless fermions considered 
in this paper. This model has been introduced by a number 
of researchers in different forms Consider a 

Hamiltonian for mass m ID spinless fermions. 
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with a zero-range interaction potential Vf{z) that has a non- 
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trivial action on the space of antisymmetric wavefunctions 
\\tf{xi, . . . ,xn). The zero-range interaction potential Vf{z), 
parameterized by an interaction strength ^{q, can be chosen 
in a number of ways, including | ^ 



and 12 



-0 2^ \pgw ^ 



(3) 



At first glance, it is not obvious that these and the other |7, 
|5,'To', 28] zero-range potentials for ID spinless fermions are 
equivalent. However, when treated exactly, they all lead to 
the same logarithmic derivative of the relative wavefunction 
where two particles touch (z is the distance between 
the particles) |0]: 



d\\freliz) 1 
dz \\trel{z) 



(4) 



The dimensionless interaction strength for this model is jf = 
'"^2°" SO that weak interactions (jf ^1) occur for low density 
and small and strong interactions (jf ^1) occur at high 
density and large ^(q. 

Recently, one of us (B. Granger), in collaboration with 
D. Blume, showed that this model of ID fermions with 
zero-range interactions describes the physics of spin-polarized 
fermionic atoms in quasi- ID harmonic traps at ultracold tem- 
peratures |8]. Then, the effective ID interaction strength 
depends on the ratio Vp/a^j_, where Vp is the p-wave scatter- 
ing volume of the atoms and a± the length scale of the tight 
harmonic confinement. Thus, by changing either Vp (using a 
Feshbach resonance |29]) or a± (by modifying the trap) 
can be tuned over a large range of values, including infinity. 
Thus, it should be possible to study strongly interacting ID 
fermions using ultracold atoms; just as strongly interacting 
ID bosons have been studied in recent experiments 1 14, 15]. 

To calculate the single particle correlations in the strongly 
interacting limit {jf °o) we map the ID fermions to a sys- 
tem of ID bosons with reciprocal interaction strength. This 
mapping was first used by Girardeau |4| to relate impenetra- 
ble bosons to noninteracting fermions. Recently, the map- 
ping has been extended to relate ID fermions with any in- 
teraction strength to ID bosons with 5-function interactions 
More specifically, any antisymmetric wave- 
function , . . . ,ZAf) for ID fermions with zero-range 
interactions of strength is related to a symmetric wave- 
function , . . . ,zn) of ID bosons with 5-function inter- 
actions of strength through the mapping |j4|] : 



\|/f(zi,...,za?) 
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When the interaction strengths of the bosons and fermions are 
chosen to be reciprocals of each other. 
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the energy spectrum, thermodynamic properties and local cor- 
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Yl'jZk ^g^izj — Zk) = 1) of the two systems are 
identical. Thus, for example, strongly interacting fermions 
(y^ oo) share these properties with the noninteracting ID 
Bose gas. 

Not all properties, however, are shared by the fermions and 
their bosonic partners. In fact, any nonlocal correlation func- 
tion will be different. Thus, nonlocal correlations provide an 
important and unique signature of the actual exchange statis- 
tics of the physical particles. We now derive the single particle 
correlations and momentum distributions for the ID fermions 
described above in the strongly interacting limit (y/r — > oo). 

General form of the density matrix. For both bosons and 
fermions, the single particle reduced density matrix is defined 
as 

rh 

P(Z,Z') = A^ / y\^{z,Z2,---,ZNW{z!,Z2,...,ZN) 
J a 

X dZ2---dZN, (7) 

SO that Tr [p] — N. The ground state wavefunction, \\tB, for A^ 
noninteracting bosons has the familiar form 



^b{zu---,zn) 



-flUzi 



(8) 
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where (|)o(z,) is the normalized ground state solution to the 
single particle Schrodinger equation in the region a < z < b. 
From the form of \\tB in Eq. (|8}, the reduced density matrix 
for noninteracting bosons, Pb{z,z'), is easily found and quite 
simple: 



PB{z,z')=NUz)roiz'). 



(9) 



Finding the density matrix, Pf for ID fermions with 
strong (jf oo) zero-range interactions, on the other hand, 
is nontrivial, but can be accomplished by using duality. We 
construct the interacting fermion wavefunction in the jf ^ ao 
limit using the noninteracting boson wavefunction, Eq. (|8j, 
and the mapping, Eq. (|5}. Then, substitution of the fermion 
wave function \\tf, Eq. into the definition of the density 
matrix, Eq. Q, yields an integral expression for the single 
particle fermionic density matrix in terms of the orbitals of 
the noninteracting bosons: 

Pf(z,z') = NUzWoiz') J' (jlsgniz-z,)sgn{z' ~z,)j 

X n (^8ni^j-^k)f^(^i^kWoi^k)\ 
\l<i<k / 

X dz2---dzN- (10) 

Noticing that {sgn{zj — Zk))^ = 1 in this expression, the multi- 
dimensional integrals can be separated into a product of A^ — 1 
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integrals over each dzi- A careful treatment of these inte- 
grals gives the final expression for the single particle reduced 
density matrix for ID fermions in the strongly interacting 
(y^ ao) limit: 



PF{z,z') = PBiz,z'){l-2P{z,z')) 



In this expression 



N 



PB{z",z")dz" 



(11) 



(12) 



is the probability of finding a noninteracting boson between 
z and z'. This result gives a completely general relationship 
between the single particle density matrices of noninteracting 
bosons and strongly interacting fermions in ID. 

In many cases it is possible to take the thermodynamic limit 
of these results (N °°, L ^ °o but n — N/L fixed, where 
L is some characteristic length for the system). If P{z,z!) = 
LP{z,z!) is independent of and L, then the expression 



: lim (1 

N^oa 



x/Nf 



(13) 



can be used to take the thermodynamic limit of the fermion 
density matrix, Eq. ( II 1> : 



-2nP(z,z!] 



(14) 



Typically P{z,z') ^ \z — z!\, so that the ID fermions gener- 
ally exhibit an exponential decay of correlations with distance 
compared to noninteracting bosons. 

Periodic system. Now that we have the general expression, 
Eqs. ( II 1> and M2\ . for the single particle density matrix in 
the limit y/r oo, we can apply this result to specific exter- 
nal potentials. First, consider a constant potential with peri- 
odic boundary conditions over the interval [0,L]. The ground 

state solution for each orbital is just (j)o(zj) = \J\, so from 
Eq. (H, Pb = n. The fermion density matrix then follows from 
Eqs. dlD and dEJ: 



Pf(z,z')-«(l-2 

The thermodynamic limit (A^ - 
fixed) of this result. 



(15) 



L oo and n = N/L 



pf (z,z') =ne 



-2nh'-z\ 



(16) 



exhibits the promised exponential decay with distance. Al- 
ternatively, the single particle correlations can be charac- 
terized by the correlation function, gi{z) — Pf(z,0)/« = 
exp(— 2n|z|), which is normalized as gi{0) = 1. This expo- 
nentially decaying single particle correlation function is one 
of the central results of this paper and clearly shows the strik- 
ing effect that the boson^fermion mapping, Eq. (|5}, has on 
the nature of single particle correlations. More specifically, 
the long-range order of T = 0, noninteracting ID bosons be- 
comes an exponential decay of correlations (with correlation 



length ^ = ^) in the strongly interacting ID fermions con- 
sidered here. Furthermore, this behavior is a dramatic depar- 
ture from the power law decay of correlations exhibited by the 
r = 5-interacting ID Bose gas and Luttinger liquid. 

Of particular importance for experiments with ultracold 
atoms is the momentum distribution n{k), which is obtained 
as the Fourier transform of the correlation function gi{z) : 



n{k) 



1 

27t 



e"'^giF{z)dz 



2n 



K{4-n^ + k^) ' 



(17) 
(18) 



This Lorentzian momentum distribution shows the signature 
effect of strong interactions in a system of zero-temperature 
fermions. That is, interactions have completely smoothed out 
the sharp decrease of the momentum distribution at the Fermi 
momentum kp w Inn. 

It is also interesting to note that the width of this momen- 
tum distribution is given by the linear density n. As « ^ this 
momentum distribution becomes identical to that of a nonin- 
teracting gas of zero-temperature ID bosons, b{k). This is 
quite an unexpected momentum distribution for a strongly in- 
teracting gas of fermions. While it is tempting to speculate 
about the physical meaning of this feature (is it a noninteract- 
ing Bose gas of molecules?), we warn that our results are only 
relevant in the high density limit (jp oo) limit. 

Harmonic confinement. Finally, we can determine the den- 
sity matrix and momentum distribution for ID fermions of 
mass m confined to a harmonic oscillator of frequency CO^ in 
the z direction. The normalized orbitals are 



Mzj) 



1 



(19) 



where a^ = y/h/moi^. The density matrix for noninteracting 
ID bosons, Eq. (|9}, is then 



Pb = ne 



(20) 



where we define the linear density for the harmonic oscillator 
as « = ^r^S ■ density matrix for our strongly interacting 
ID fermions, Eqs. dl 1> and (I12> . can be written in terms of 

Erf{x)^j^S^e-''dt: 



PHz,z')=Pb(1- I £'•/(-)-£'-/(-) I) 



N-l 



(21) 



From this expression, we calculate the momentum distribution 
n{k) numerically, by Fourier transforming the density matrix. 



PF{k,k') 



1 
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e'kze-''''''p{z,z')dzdz', (22) 



so that n{k) = Pf{k,k). 

Figure ^ shows the calculated momentum distributions of 
harmonically confined, strongly interacting ID fermions for 

= 5, 20 and 100 particles. The main figure shows clearly 
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interactions. 
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Figure 1: Momentum distributions n{k) are shown as a function of k 
for a ID gas of strongly interacting spinless fermions with ID har- 
monic confinement of frequency to,. The main figure shows the 
momentum distributions, normalized to unity with the momentum 
given in units of I /a,. The inset shows the momentum distribution 
(normalized to N) with the momentum given in units of the density 
n = N /a^y/n. For the N = <x> case we have plotted the Lorentzian, 
Eq. <18> . that corresponds to the thermodynamic limit. 



that the widths of these distributions increase in proportion 
with the Hnear density n. This is the opposite behavior from 
what has been calculated for harmonically confined ID inter- 
acting bosons i25ll . The oscillations at small momentum are 
due to the harmonic confinement rather than the details of the 



At a more quantitative level, the large momentum tails of 
these distributions can be understood by taking the thermo- 
dynamic limit of the density matrix, Eq. (12 1> . This limit is 
specified by letting a^-y/jt °° and N ^ °° with n = N/a^^/n 
fixed. In this limit, we can use the small argument ap- 
proximation Erf{j^) « and Eq. il3\ to obtain p/r = 

PBexp(— 2n|z — z'l). Thus, in the thermodynamic limit, the 
short range correlations in a harmonic potential are identical 
to those of a periodic system. In the momentum distributions, 
this shows up as a Lorentzian tail, Eq. jl8L with a width of 
4«, where n is the linear density. This feature is confirmed in 
the inset of Fig.^ where the momentum is plotted on a log- 
arithmic scale in units of the density n to show how at large 
momentum all the distributions line up with a Lorentzian of 
width 4« (solid line). This Lorentzian momentum tail will 
be a clear experimental signature of the strongly interacting 
regime discussed here. 

In conclusion, duality in ID many body systems provides 
a powerful method for understanding and calculating corre- 
lations in the strongly interacting limit. The system we con- 
sider, namely, strongly interacting spin-polarized ID fermions 
with zero-range interactions, is dual to a noninteracting gas of 
bosons and provides an excellent prospect to study duality in 
actual experiments with ultracold atom gases. The exponen- 
tially decaying correlations described here will have a clear 
and measurable effect in these experiments and will hopefully 
spur on further theoretical inquiries as well. 

This work was funded by the NSF under grant PHY- 
0354882. 



[1] S. Tomonaga, Prog. Theor Phys. 5, 544 (1950). 

[2] J. M. Luttinger, J. Math Phys. 4, 1 154 (1963). 

[3] L. Tonks, Phys. Rev. 50, 955 (1936). 

[4] M. D. Girardeau, J. Math Phys. 1, 516 (1960). 

[5] S. Coleman, Phys. Rev. D 11, 2088 (1975). 

[6] T. Cheon and T. Shigehara, Phys. Rev Lett. 82, 2536 (1999). 

[7] S. Albeverio, L. Dabrobski, and S. M. Fei, 

arXiv quant-ph/0001089 (2000). 
[8] B. E. Granger and D. Blume, Phys. Rev. Lett. 92, 133202 

(2004). 

[9] M. D. Girardeau and M. Olshanii, arXiv cond-mat/03093961 

(2003) . 

[10] H. Grosse, E. Langmann, and C. Paufler, 

arXiv math-ph/0401003 (2004). 
[11] A. Goiiitz et al., Phys. Rev. Lett. 87, 130402 (2001). 
[12] S. Dettmer et al., Phys. Rev. Lett. 87, 160406 (2001). 
[13] S. Richard etai, Phys. Rev Lett. 91, 010405 (2003). 
[14] B. Paredes et al.. Nature 429, 277 (2004). 
[15] T. Kinoshita, T. Wenger, and D. S. Weiss, Science 305, 1125 

(2004) . 



[16] M. Olshanii, Phys. Rev. Lett. 81, 938 (1998). 
[17] E. H. Lieb and W. Liniger, Phys. Rev. 130, 1605 (1963). 
[18] J. B. McGuire, J. Math Phys. 5, 622 (1964). 
[19] C. N. Yang and C. P Yang, J. Math Phys. 10, 1 1 15 (1969). 
[20] S. Inouye etal.. Nature 392, 151 (1998). 
[21] H. G. Vaidya and C. A. Tracy, Phys. Rev Lett. 42, 3 (1979). 
[22] D. M. Gangardt and G. V. Shlyapnikov, Phys. Rev. Lett. 90, 
010401 (2003). 

[23] K. Y Kheruntsyan, D. M. Gangardt, P. D. Drummond, and G. V. 

Shlyapnikov, Phys. Rev. Lett. 91, 040403 (2003). 
[24] M. Olshanii and V. Dunjko, Phys. Rev. Lett. 91, 090401 (2003). 
[25] G. E. Astrakharchik and S. Giorgini, Phys. Rev A 68, 031602 

(2003). 

[26] F D. M. Haldane, J. Phys. C 14, 2585 (1981). 

[27] K. Kanjilal and D. Blume, arXiv:physics/0406022 (2004). 

[28] M. D. Girardeau, H. Nguyen, and M. Olshanii, 

arXiv cond-mat/0403721 (2004). 
[29] C. A. Regal, C. Ticknor, J. L. Bohn, and D. S. Jin, Phys. Rev. 

Lett. 90, 053201 (2003). 



